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Abstract 

We consider a driven diffusive system with two types of particles, A and B, coupled at the 
ends to reservoirs with fixed particle densities. To define stochastic dynamics that correspond 
to boundary reservoirs we introduce projection measures. The stationary state is shown to be 
approached dynamically through an infinite reflection of shocks from the boundaries. We argue that 
spontaneous symmetry breaking observed in similar systems is due to placing effective impurities 
at the boundaries and therefore does not occur in our system. Monte-Carlo simulations confirm 
our results. 
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I. INTRODUCTION 



Systems of driven diffusing particles attract attention because, despite their relative sim- 
plicity, they embrace a whole range of critical phenomena far from thermal equilibrium 
[1-3]. One of the remarkable features is the appearance of phase transitions induced by 
spatial boundaries of a system, studied in detail for models with one species of particles 
[4-8] and for some multi-species models [9-12]. The ability of a nonequilibrium system to 
"feel" the boundaries constitutes a key feature of driven systems: in fact, the boundaries 
dominate the bulk giving rise to the phase transitions. Therefore it is clear that boundary 
conditions play a crucial role and have to be chosen with care, strictly adequate to the 
physical situation which is being modelled. To this end one has to define specific boundary 
rates for injection and extraction of particles at the ends of the system. There are several 
possible strategies. 

One of them is to postulate the most simple boundary rates, and treat the rates themselves 
as parameters. A system is studied, then, as a function of those parameters, so that one 
tries to keep their number small. This approach is used, e.g., in [4, 9]. 

Another approach, let us call it particle boundary reservoir approach, treats the boundary 
problem of a driven diffusive system as if it was coupled at the ends to reservoirs of particles 
with fixed particle densities. In this case, the rates are dictated by the reservoir particle 
densities and by themselves do not, in general, look simple. The boundary rates are fixed 
such that the stationary measure of the process with equal reservoir densities on the left 
and on the right is essentially equal to the stationary measure of the infinite system with 
the same particle densities, see below for a detailed discussion. Thus reservoir particle 
densities are parameters, namely, reservoir densities of each species of particles at each 
boundary. Correspondingly, in the one-dimensional case, where there are two ( left and 
right) boundaries, the number of parameters is equal to 2*number of species in the system. 
This approach is used in [8, 10, 13]. We argue that in this case there exists a well-defined 
hydrodynamic limit, and the original stochastic problem has correspondingly a well-behaved 
coarse-grained description via conservation law equations. The latter constitutes a powerful 
tool to study the physics of the system, e.g., phase diagrams and the location and type of 
phase transitions. 

In some important cases the two approaches are equivalent. For instance, for the ASEP 
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with open boundaries, the rates of injection a and extraction (3 respectively are expressed as 
a = pl, (3 = 1—pR in terms of the particle boundary reservoir approach. In more complicated 
settings one may attribute effective boundary densities to given boundary rates. In this way 
a very large class of single-species systems with open boundaries can be described and their 
phase diagram can be predicted in terms of the boundary densities [8] . For attractive single- 
species systems the validity of this approach has been proved rigorously in very recent work 
[14]. The purpose of this paper is to demonstrate that in two-species models, however, the 
two approaches cannot generally be substituted with one another. We exhibit this difference 
and discuss its significance for the appearance of spontaneous symmetry breaking in an ABO 
model [9], which is a paradigmatic two-species model [15]. 

To this end, we first introduce the model (Sec. II), discuss its hydrodynamic properties 
under Eulerian scaling (Sec. Ill) and show how to construct appropriate boundary reservoirs 
(Sec. IV). The strategy employed here can be applied to any stochastic particle system. 
Then we show (Sec. V) that the phenomenon of infinite reflections of shock fronts from 
the boundaries [16] exists also in the ABO system which is much studied, and is known to 
have many intriguing features. Using this we arrive at the conclusion (Sec. VI) that the 
so-called bridge model [9] (which is an ABO model with specific injection and extraction 
rates) has spontaneous symmetry breaking due to the fact that the boundary injection and 
extraction rates place effective obstacles at the ends of the chain, rather than corresponding 
to boundary reservoirs with some effective boundary density. We demonstrate that with 
these effective obstacles removed, there is no symmetry breaking phase transition. Thus, 
symmetry breaking can be viewed as an impurity- induced effect. Moreover, we argue that the 
original bridge model may have symmetry breaking transition precisely at the point where 
a hydrodynamic description with effective boundary densities breaks down. We summarize 
our main conclusion in the final Sec. VII. 

II. THE MODEL 

The ABO model is defined on a chain of length N. There are particles of type A and of 
type B, distributed along the chain. Any site of the chain can be empty, or contain at most 
one particle. In the bulk, exchange with the nearest neighbours happens, e.g., a particle A 
can exchange with a hole on its right AO =>- OA with rate T^o, etc.. The complete set of 
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rates is given as follows 



AO =>• OA 


with rate T^o 


r\ A \ A r\ 

OA =>- AO 


with rate 1 a 


BO OB 


with rate T^o 


OB =>■ SO 


with rate T b 


AB =>■ £A 


with rate Tab 




with rate T Bj 4 



(1) 



The boundary conditions correspond to reservoirs of particles with fixed densities pf,p£ 
and //ft, at the left and at the right end of the chain respectively and are defined in the 
section IV. For a specific choice of rates [17], 

Tab — ^ba = (Tao — ^oa) + (Tob — ^bo), (2) 

the stationary state on a ring is given by a product measure, which does not crucially affect 
the qualitative properties of the system, but significantly reduces computational difficulties 
(since no approximation in the calculation of the bulk flux, boundary conditions are re- 
quired). We shall follow this choice below, making comments on the general situation when 
necessary. In most cases, we shall take 

Tao = Tob = Tab/2 = 1, T a = T B o = T B a = (3) 
without losing generality. 

III. HYDRODYNAMIC LIMIT EQUATIONS 

Here we shall treat for simplicity the case (3) where there is a product measure state on 
a ring. From this property and the hopping rules (3) exact stationary fluxes are given by 

Ap a ,p b )=p a (i-p a + p b ) (4) 

f{p\p B ) = -p B {l-p B + p A ). (5) 

Following [18], the hydrodynamic equations of our process may be obtained by averaging, 
factorizing and Taylor expanding the lattice equations of motion with respect to the lattice 
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constant a< 1 . In the first order of the expansion, one obtains 



dp A dj A (p A ,p B ) 



d_ 

dx 




df dx 



dp* , df(p A ,p B ) 



d_ 

dx 
d „ 



(6) 



dt' dx 




For a finite system with open boundaries this is supplemented with the boundary conditions 



The system of equations (6), after the substitutions p = 1 — p A — p B , u = p B — p A , becomes 



which in its inviscid form (e = 0) has been derived rigorously for our model defined on the 
infinite lattice [19] and is known in the theory of conservation law equations as the Leroux's 
equation. It appears in many applications, further details about this equation can be found 
in [19, 20]. Note that the viscosity term is required for uniqueness of the solution of the 
hydrodynamic equation. For an infinite system the choice of viscosity has some arbitraryness. 
However, in the presence of the boundaries the exact form of the viscosity matrix becomes 
important in systems with more than one conservation law [21]. The question how the 
macroscopic boundary condition (7) on the PDE can be realized on microscopic lattice gas 
scale is addressed in the following sections. 

IV. BOUNDARY RESERVOIRS AND PROJECTION MEASURES 

In driven (nonequilibrium) systems, the steady state is governed not only by local inter- 
actions in the bulk, but also by interaction with boundaries. To define a boundary reservoir, 
we introduce a projection measure for a subsystem of a larger system (defined on a spatial 
domain fi) by defining the following marginal: a measure is a projection measure for a 
spatial domain Q C Q if all its correlations are identical to those of the whole system Q 
inside the domain Q. 



p z (0,t)=pl, p z (N,t)=p z R , Z = A,B 



(7) 




(8) 
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In what follows we consider as domain Q a finite system of N sites as a subsystem 
of an infinite system where Q is the integer lattice Z. For illustration we first consider 
N = oo (semi-infinite system with a left boundary). In the case of a product measure, the 
projection measure is especially simple: it remains a product measure for all sites in Q, i.e., 
homogeneous initial distribution with an average densities Pb u ik> Pbulk °^ ^ anc ^ ^ particles. 
In vector notation [1] the projection measure for the semi-infinite system is given by 



\m) = n 



k=l 



Pbulk 
Pbulk 



(9) 



1 _ n B - n A 
\ L Pbulk Pbulk/ 

Similarly, for a finite subsystem of the length iV one has the projection measure 

/ \ 

N 

ip>=n® 



i=i 



u 

V 

1 — u — V 



(10) 



Pbulk' V — Pbulk- 



where u 

So far we have discussed only measures. To make the connection to stochastic dynamics 
we assume that the particle system defined on the infinite lattice has a product measure 
as its steady state. In order to define dynamics for an open system with creation and 
annihilation of particles corresponding to a boundary reservoir with densities pf = u, pf = v 
for a semi-infinite system, we demand that the projection measure (9) stays invariant for 
u = Pb ulk , v = Pbuik- This is a natural requirement, since one has a perfect match of the 
product state with the boundaries. In the same way, we can define the right boundary 
reservoir for a semi-infinite chain k = — oo, 0, and for a finite chain of the length N, with 
the projection measure (10), where pf, = Pr = u i Pl = Pr = v - Some examples illustrating 
more general situations (e.g., involving an Ising measure) can be found in [8, 10, 13]. 

In general, there are many ways to obtain boundary rates which leave the projection 
measure invariant. In order to calculate boundary rates for our model we first require that 
particles are injected only at the edges (sites 1,N) of the lattice, corresponding to the 
nearest neighbor hopping in the bulk of the chain. Furthermore we note that, in particular, 
one requires the stationary flux in the system with equal reservoir densities Pl = Pr = M > 



Pl = Px = v to be equal to the stationary flux in an infinite system with the densities u, v. 
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The stationary flux is 

j A (u,v) = u(l -u + v) (11) 
j B (u,v) = -v(l -v + u). (12) 

The different signs of j A and j B are due to the fact that A-particles and 5-particles hop in 
opposite directions. On the other hand, the stationary flux of particles A through the left 
boundary is 

j£ = (0)T L A0 + (B)T L AB (13) 

where Y^ Y is the rate of exchange of particle of the sort X (from the left reservoir) with 
the particle of sort Y at the first site and (Y) is the stationary expectation value to find a 
particle Y at the first site. In the stationary state, the latter is equal to the density (see 
(10)): (B) = v , (0) = 1 — u — v. Inserting these expressions in (13), we get 

3$ = {l-u-v)T L AQ + vT L AB 

Repeating the arguments for the flux of particles B through the left boundary, we obtain 

jE = — (B)Y L B - (B)T L AB = -v{Y L 0B + Y L AB ) 

We require the above expressions to be equal to the stationary fluxes (11), (12). This leaves 
one parameter undefined, since there are 3 unknown boundary rates Y and two relations. To 
get an additional relation, consider another correlation, e.g. the probability to find particles 
A, B on the first two sites, (AB). The time evolution of this correlation is 

|-(AB> = Y L A0 (0B) + Y 0B (A0B) + Y AB (AAB) - Y AB (AB). 
ot 

In the stationary state, left hand side vanishes, and the correlations are given by (10) 
(OB) = (1 - u - v)v, (AB) = uv, (AAB) = u 2 v, (AOB) = u(l - u - v)v. Substituting 
the bulk rates Y X y from (3), we get = u, which, together with the relations j A = j A , 
3l — J B fi xes an boundary rates to be 

Y L A0 (u,v) = u; Y L AB (u,v) = 2u- Y% B (u, v) = (1 - u - v). (14) 

It is interesting to note that the rates (14) can be obtained from another simple argument. 
In a boundary reservoir with the densities p A = u, pf = v, the probability to find a particle 
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A, B and a hole is p(A) = u, p(B) = v, p(0) = l — u — v. We assume the rates of hopping 
of the particles from the reservoir to the site 1 of the system to be given by the bulk rates 
(3) multiplied by corresponding probabilities, i.e., 

T L A0 = p(A)T AQ , r£ B =p(0)r„B, T L AB =p{A)T AB . (15) 

This yields (14). 

Analogously, on the right boundary, the rates of exchange of particle X at site N 
with the particle Y (from the right reservoir) is 

T%=p(0)T AO , T^ B =p(B)T 0B , T%=p(B)T AB (16) 

or 

r%(u,v) = (l-u-v), r* B (u,v)=v, T%{u,v)=2v. (17) 

One verifies straightforwardly that for the choice of the rates (14), (17) the projection 
measure (10) is stationary. 

For the boundary reservoirs with arbitrary densities p^, p^ on the right and pf, pf on 
the left we shall use the boundary rates given by r§ y (p^, p^) and T^ y {Pl-> Pl) respectively 
We shall call such boundary rates PDE-friendly, since, as we show in the next section, they 
allow for an adequate hydrodynamic description for a finite open system. 

V. INFINITE REFLECTIONS OF SHOCK FRONTS FROM THE BOUNDARIES. 
ABSENCE OF SPONTANEOUS SYMMETRY BREAKING. 

It is an open problem how to obtain a hydrodynamic description of a two species system 
with open boundaries. A step towards such a description has been made in [16]. It has been 
demonstrated that the evolution to the stationary state can proceed through an infinite 
sequence of reflections of a domain wall with the boundaries of the system. It is interesting 
to study the reflections for the system (3) to see if it provides another example where this 
scenario applies. Note that the model we investigate here and the one investigated in [16] 
have topological differences: in the latter, particle dynamics on two parallel chains was 
considered, and particles on both chains hopped in the same direction. Instead, in the 
presently considered model (3), particles occupy one chain and hop in opposite directions. 
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A. Reflection maps 



To pose a reflection problem, take a half-infinite chain with only one, e.g., left boundary, 
corresponding to left boundary reservoir densities p A = u, p£ = v , and choose an initial state 
to be a projection measure of a stationary state for an infinite system onto the positive half- 
infinite chain (9), i.e., homogeneous initial distribution with an average densities Pb u ik> Pbulk 
of A and B particles. We consider bulk rates (3) and corresponding boundary rates (14). 

If the bulk densities coincide with the left boundary densities p^ ulk = u, p B u y k = v, the 
initial distribution (9) stays invariant because of the perfect match of the initial state with 
the boundaries. Otherwise, there is a mismatch at the boundary which has to be resolved. In 
this case, Monte Carlo simulations lead to one of the following scenarios: (a) a thin boundary 

A B 

layer develops interpolating between the bulk densities p bulk and boundary densities u, v, and 
stays always attached to the left boundary, (b) a shock wave of the densities r A , r B develops 
and propagates to the right, (c) a rarefaction wave with the density r A , r B forms and spreads 
to the right. The scenario (b) is demonstrated on Fig.l, showing Monte-Carlo evolution of 
a state (9) with p^ = (0.1, 0.7), u, v = 0.25, 0.55 after t = 300 Monte Carlo steps. Like in 
[16], we mimic the half-infinite chain with a finite chain with the right boundary conditions 
matching perfectly the bulk densities. As a result, after time t the density profile, while 
staying unchanged on the right boundary (see Fig.l), at the left boundary develops the 
shock with the densities r A m 0.348, r B p& 0.131 spreading to the right. Due to particle 
number conservation in the bulk, a Z-component of the shock interface moves with the 
velocity 

yZ = Jbulk Jr (lg) 
Pbulk - r 

where we used the short notations j£ uJk = j z (p A u \ k , Pbulk) > Jr = j z (r A ,r B ). Since there is 
an interaction between A and B species, the velocities must coincide in both components: 
V A = V B , or 

A A „*.A ~B ~B 

Jbulk Jr _ ./bulk Jr (10} 
n A _ r A~ n B _ r B' V iy J 

Pbulk ' Pbulk ' 

defining implicitly the allowed location of the points r A ,r B . In our case, the above relation 
becomes (see 4,5): 

Pb A uik(l " Pb^uik + PU) -r A (l-r A + r B ) _ -pf uIk (l - p* Ik + p A nlk ) + r B (l - r B + r A ) 



bulk ' Pbulk ' 



(20) 
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Here we implicitly assume the fluxes of the reflected wave to be given by stationary fluxes 
with the densities r A ,r B . This supposition is justified because for random initial conditions 
used here the system away from the shock is locally stationary. 

Analogously one can consider the right boundary reflection, taking a half-infinite chain 
with right boundary only and initial product measure state (9) where the product goes from 
— oo to L. An example of a right boundary reflection is shown on Fig.2. Again, the shock 
wave resolving the mismatch at the boundary appears and spreads to the left. For the same 
reasons as were given above for left boundary reflection, the densities of reflected waves 
r A ,r B satisfy relation (20). Solutions of (20) are two families of straight lines 

r B -Pbuik = (r A -p^i k )7 + (Pbui k , Pbuik), 7+(s, *) = ^ (-2 + s + t + y/(s + t -2)* - 4st) . 

(21) 

^-Pbuik = (r A -Pbuik)7-(Pbuik, Pbuik), 7-(s, t) = y s (-2 + s + t - v /( s + t-2) 2 -4st) . 

(22) 

The solutions with bigger 7 corresponds to the left reflection, and another to the right 
reflection. Note that 7±(s, £) are real in the physical domain < s + 1 < 1. The meaning 
of the above equations is twofold. One of them was given above, that is, for any given 
Pbuik' Pbuik' the reflected waves must have densities satisfying (21), (22). Another, alternative, 
is: given the reflected wave densities r B and r A , the initial bulk densities ( which have led 
to those reflected waves), must satisfy the equation equivalent to (21,22). 

pU-r B = ( P U-r A ) 1± (r A ,r B ), (23) 

where 7±(s, t) are defined in (21,22). The equation above is a solution of (20), because 
the latter is completely symmetric with respect to the exchange r z <-> p^ ulk . On the ba- 
sis of Monte-Carlo calculations and investigation of the hydrodynamic limit equations the 
following observation can be made: 

The collection of all possible densities r A ,r B of the waves, resulting from the left boundary 
reflection (with boundary densities u,v) with whatever initial conditions Pbuik>Ptmik> consti- 
tute some curve C(u,v), completely parametrized by this point u,v, and containing this 
point. All initial bulk densities satisfying (23) for certain point r A ,r B of C, result in the 
reflected wave with the densities r A ,r B . The same is true for right boundary reflection, with 
the curve C replaced by TZ. 
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The mapping from arbitrary p^ ulk , Pbuik to resulting r A , r B is described by means of reflec- 
tion maps introduced in [16]. Typical examples of left and right reflection maps are shown 
on Figs. 4 and 3. 

Each of the reflection maps of the type shown in Figs. 3 and 4 is parametrized by 
a single point u, v (which has a meaning of the left/right boundary density for left/right 
map). Each point < u + v < 1 corresponds to a different map. Classification of all 
different maps requires separate investigation. Here we shall find out what happens if we 
fix the right and the left boundary corresponding to the reflection maps Fig. 4 and Fig. 3. 
Let us take an initial homogeneous state with the densities, fitting the density of the left 
reservoir, p^ u i k , p^ u \ k = Pl^Pl- As the result of the first interaction with the right boundary, 
a reflected wave forms, with the densities 7q , Vq given by the intersection of the line (22) 
with the curve 1Z of Fig. 4. This reflected wave hits the left boundary then, resulting in the 
next reflected wave with the densities rj 4 ,rf . The left reflection is controlled by the curve C 
at Fig. 3, therefore the densities of the reflected wave rj 4 , rf are given by the intersection of 
the line (21), where new p^ ulk , p^ ulk = rjf , rjf , with the curve C. Since the curves 1Z,£ do not 
coincide with any of the lines defined by (14), (17), this process continues forever, though 
converging to the stationary state S, the point of the intersection of C with 7Z, see Fig. 5. 
After 2k — 1 reflections, the wave with the densities r^_i, r 2k-i wm hit the right boundary, 
producing the reflected wave of the densities r^., r^., corresponding to the intersection of 
the line (22) ( with Pbuik> Pbuik = r 4-i> r 2k-\) with the curve 1Z. The latter wave, hitting 
the left boundary, produces the next reflected wave, and so on. This process is depicted 
schematically on Fig. 6. Several remarks are in order. 

• The densities {r£, rf }^L of the reflected waves constitute a converging sequence. 
Convergence is exponential, 5u n+ 2k — e~ Kh 8u n , for n — > oo where we denote by 5u n 
the deviation from the stationary density at the n-th step. e~ K < 1 is a constant 
depending on tangential derivatives a, (3, a±, (5\ of the curves 1Z,£ and the character- 
istic curves(14), (17), respectively, at the stationary point S, e~ K = (1 — )(1 — 

iB)^^ 1 + t ~$m" > ( 1 ~ WW ))■ Note that if at least one of the characteristic deriva- 
tives happens to coincide with «i, fli(a = a.\ or (3 = Pi) then the sequence converges 

in one step. 

• The velocities of reflected waves for large k converge to the finite characteristic ve- 
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locities computed at the stationary point S, as eigenvalues of the flux Jacobian, see 
[18]. 

• The stationary densities are not reached at any finite time. Strictly speaking, the sta- 
tionary density is reached exponentially in time with characteristic time proportional 
to the length of the system N. On the contrary, in one-species models like the ASEP, 
the corresponding characteristic time is of order of l/(rate of hopping). 

We conclude that for the system (3) with PDE-friendly boundary rates 
^xy(Pli Pl)i ^xy(Pri Pr)i given by (14), (17), the convergence to a stationary state 
through an infinite number of reflections is a generic feature of the dynamics. This is 
confirmed by the good agreement of the Monte-Carlo simulation of the lattice model with 
the numerical integration of hydrodynamic PDE's (6), see Figs.l, 2, 5, 7. Without the 
PDE-friendly boundary rates already the result of the first reflection would be unpredictable. 



B. Symmetric boundary conditions 

In the preceeding subsection we gave an example for the time evolution of the system in 
case of generic boundary conditions. Here we consider a situation when the boundary rates 
possess the symmetry with respect to a simultaneous right-left reflection and an exchange 
of A particles with B particles, i.e., the same symmetry as in the similar models where 
spontaneous symmetry breaking has been observed [9, 22]. This amounts to the following 
restriction of the boundary reservoir densities: 

Pr = Pl> Pr=Pl- ( 24 ) 

Analysis of reflection maps for our PDE-friendly boundary rates (14,17) allows us to con- 
clude, that there is no spontaneous symmetry breaking — the stationary state is always 
symmetric (i.e., corresponds to equal densities of A- and 5-particles). A typical reflection 
map for the case (24) is presented in Fig. 7. The reason why there is spontaneous symmetry 
breaking in the bridge model [9] will be given in the next section. 
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VI. WHY THE HYDRODYNAMIC LIMIT FAILS FOR THE BRIDGE MODEL 



The bridge model is defined on a finite chain of the length N, where the processes in the 
bulk with rates 

= T 0B = 1, Tab — Q, ^oa = Tbo = ^ ba = (25) 
are complemented with the processes at the boundaries. On the left boundary site 

r io = a = aT Ao, Tq B = (3 = (3T 0B , T L AB = (26) 

and on the right boundary site, 

r* B = a = ar 0B , T% = P = f3T A0 , V% = 0. (27) 

The model is invariant with respect to simultaneous left-right and A -B interchange. Nev- 
ertheless, for certain range of parameters, typically for large a and small /3, the state of the 
system is characterized by a phase where the average densities of particles are not symmet- 
ric: the symmetry between A and B is broken spontaneously, and it takes a system a very 
long time (which grows exponentially with the system size) to get from the state with a 
prevalence of A-particles to the state with a prevalence of -B-particles. The exchange rate q 
is not crucial for the occurence of this intriguing phenomenon. 

(a) Let us consider q = 2 and try to see if one can associate a reservoir density as described 
above to the rates (26), (27). In terms of the probabilities p(A) , p(B) , p(Q) to find a particle 
of sort A,B and respectively in the reservoir, the left reservoir rates are given in (15). 
Comparing (26) with (15), we obtain relations 

p(A) = a,p(p) = P,p(A) = 0, 

incompatible between themselves. 

(b) Let us nevertheless try to find effective boundary densities for the rates (26, 27) for 
general q. It is instructive to analyze the case when both injection and extraction are small 
a, /? <C 1 and a < (5. Then, the typical time to inject a particle is bigger than the typical 
time to extract a particle r in = a' 1 > /3' 1 = r out , and both T in ,r out ^> 1. In this case, one 
can associate boundary densities as it is done in the usual asymmetric exclusion process 
(ASEP): — a, pf = 1 — f3, and p B = 1 — (3, p B = a. Correspondingly, the density of holes 
in the boundary reservoirs is 1 — p A — p B , yielding 

P°r = P°l=P- a. (28) 
13 



We see that for a > (3 the density of holes (28) becomes negative, which signals a breakdown 
of the reservoir picture across the point a = (3 when a,/3 <C 1. And this point lies exactly 
at the phase transition line to the symmetry broken phase [22]. Thus, we have shown that 
for parameters where symmetry breaking phase is observed, the hydrodynamic description, 
using boundary reservoirs, fails. A natural interpretation of the rates (26), (27) would be 
that they effectively correspond to impurities, placed at the left and at the right boundary 
of the system, and preventing A — B exchange. It is solely due to these obstacles, that 
spontaneous symmetry breaking actually happens. An effective hydrodynamic description 
may still be valid in the symmetric phase. 

The above considerations are valid for (a) q — 2, a, f3 arbitrary and (b) q arbitrary and 
a, (3 small. However, we believe that a similar description is valid for any q,a,(3. Since for 
q = 2 as well as for q — 1 there are no correlations between particles and vacancies, a perfect 
hydrodynamic description of the bridge model will be achieved according to the result of the 
previous section, if exchange at the boundaries between A, B happens on the same footing 
as in the bulk, see the discussion below (14), which amounts to choosing 

T L AB = V% = qa, q = l,2 (29) 

in (26), (27). For the above choice, a hydrodynamic description with fixed boundary densities 
is valid for all a, (3 leading to the complete disappearance of the symmetry broken phase 
from the phase diagram. For instance, for q = 1, the A-particle cannot distiguish between 
£>-particle and a hole, both in the bulk and in the boundaries (and the same is valid for B- 
particles), so that the system effectively separates into two one-species problem ASEP, with 
injection rate a and extraction rate a + (3. Invoking the known solution for ASEP phase 
diagram [5], [6], we find that the exact stationary state density of a PDE- friendly bridge 
model with q — 1 is 

pftat = Pftat = if a < \ 

A B 1 •( -s- 1 

Pstat = Pstat = 2' 2 

Note also that the stationary state in this case, as opposed to the q = 2 case considered 
earlier, is reached after single interaction of a shock wave, or a rarefaction wave, with the 
boundary, like in the ASEP [16]. 
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VII. CONCLUSIONS 



We have investigated the ABO model and, defining projection measures, introduced 
boundary conditions compatible with a hydrodynamic description. We have studied the 
system for boundary rates which are either symmetric or nonsymmetric with respect to 
left-right reflection and exchange of particle species. Investigating the reflection of shocks 
at the boundaries, we conclude that the system with symmetric PDE-friendly boundary 
conditions has a stationary symmetric phase in all parameter space. The relaxation to the 
stationary state proceeds by an infinite sequence of reflections with the boundaries, which 
can be described using reflection maps. Hydrodynamic limit equations are constructed and 
studied, yielding results consistent with the stochastic dynamics. For the bridge model, 
which is a special case of the ABO model, we showed that at the phase transition to the 
spontaneously broken phase the hydrodynamic description fails. While our discussion was in 
the framework of the ABO model, we believe that similar results are valid for other driven 
diffusive systems with two conservation laws, such as two-lane models [10, 18, 23, 24] or 
bricklayer models with nonconserved internal degrees of freedom [25] . 

We conclude that any choice of boundary rates which does not correspond to certain den- 
sities of boundary reservoirs is equivalent to placing an impurity at the respective boundaries. 
This entails extra complexity and rate-dependent nonuniversal behaviour. Sometimes the 
impurity thus introduced leads only to a local disturbance and correspondingly to a redefi- 
nition of boundary densities. In other cases, however, the impurity effect is highly nonlocal 
and nontrivial. This last case is very pronounced in the case of spontaneous symmetry 
breaking considered in [9]. It would be interesting to see if for any other choice of the A — B 
boundary exchange rates < T^ B = T^ B < qa there will be a region in the phase space 
a,j3 (shrinking as ^ab^ab increase) where the symmetry broken phase will exist. 
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Figure 1: Domain wall reflection from the left boundary, as given by Monte Carlo simulations. 
Particles were initially distributed randomly without correlations with constant densities p A = 
0.7, p B = 0.1, fitting the right boundary. The average profiles after 300 Monte Carlo steps are 
shown, and compared with the hydrodynamic limit (6) evolution (continuous curves). Averaging 
over 5*10 5 histories is made. The smoothing of the Monte-Carlo density profile is due to fluctuations 
in the shock position which is scaled out in the hydrodynamic description. 
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Figure 2: Domain wall reflection from the right boundary. Particles were initially distributed ran- 
domly without correlations with constant densities p A = 0.25, p B = 0.55 fitting the left boundary 
density. Average density profiles after t = 300 MCS steps are shown. The boundary densities 
are as on Fig. 1. Lines show the hydrodynamic evolution (6) for the two particle species. The 
smoothing of the Monte-Carlo density profile is due to fluctuations in the shock position which is 
scaled out in the hydrodynamic description. 
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Figure 3: Left reflection map for the left boundary densities 0.25, 0.55 of A and B particles re- 
spectively. The crossmark labelled L marks the point of the left boundary density. The curve 
containing this point L is a location of all possible densities r A ,r B of the reflected waves, except 
at the region HD. The straight lines show initial bulk densities with the same outcome of a reflec- 
tion (point of intersection of the straight line with the curve). In the HD region the initial wave 
stays glued to the left boundary. In the RW region (the upper triangle) the reflection results in a 
rarefaction wave (scaling as (x/t) with time), converging to point L. The dotted line r A + r B = 1 
shows the boundaries of the physical region. 
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Figure 4: Right reflection map for the right boundary density 0.1,0.7 of A,B particles. The 
crossmark labelled R marks the point of the right boundary density. The thick curve is a location 
of all possible densities r A ,r B of the reflected waves, except at the upper left triangle and in 
region RB. Straight lines show initial bulk densities with the same outcome of a reflection (point of 
intersection of the line with the thick curve) . In the upper corner the initial wave stays glued to the 
right boundary. In the region RB an intermediate plateau appears, with the densities corresponding 
to the broken line passing through point R. The final result of this reflection is a wave matching 
the right boundary R. The dotted line r A + r B = 1 shows the boundaries of the physical region. 
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Figure 5: Curves of the right and left boundary reflection combined together. The intersection S, 
is the stationary state density achieved through the infinite number of reflections. The crossmarks 
labelled by R and L indicate the left and right boundary densities, and the cross at the intersection 
the stationary density obtained by Monte Carlo simulation of the system of 300 sites. The system 
was equilibrated for = 4 * 10 5 Monte Carlo Steps (MCS), after which the averaging over again 
4 * 10 5 MCS and 10 different histories was done. The big star to the righthand side of S marks 
the result of a single reflection for an initial profile fitting the left boundary density point L (from 
Monte Carlo calculations). 
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Figure 6: Approach to the stationary point S through an infinite sequence of reflections. The 
filled circles show the location of subsequent domain wall densities r A , r B . The circles k, k + 2, k + 
4 correspond to the result of the left reflection, while the other circles correspond to the right 
reflection. 
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Figure 7: Curves of the right (filled triangles) and left (empty triangles) boundary reflection for 
the symmetric setting (24), as given by numerical integration of the hydrodynamic equations 
(6). Triangles labelled R and L indicate left and right boundary densities, and the point S the 
stationary density, achieved through an infinite series of reflections. The cross at the intersection 
point S marks the stationary density obtained by Monte Carlo simulation of the system of 300 
sites. The broken and dotted lines intersecting at S show characteristic curves at the point S (see 
(21),(22)). 
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